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High  performance  modeling  of  brittle  materials  is  an  efficient, 
inexpensive,  time-saving  solution  for  optimal  design  of  armor 
systems.  Quantifying  the  effects  of  high-velocity  projectiles  into 
brittle  materials  provides  an  armor  resolution  for  the  critical  need  of 
ballistic  protection  against  lethal  threats.  The  analysis  modeled  a 
cylindrical,  tungsten  carbide  blunt  projectile  into  four  confined, 
ceramic  materials  at  a  high  velocity  impact.  The  finite  element 
simulations  were  performed  using  Elastic  Plastic  Impact  Code 
(EPIC)  [Johnson  (2006)],  which  simulates  the  failure  and  particle 
breakup  of  the  target  once  the  long-rod  penetrator  strikes  at  high- 
velocity  impact.  The  history  of  the  nose  penetration  of  the 
projectile  will  be  computed  to  establish  the  most  advantageous 
design  condition  for  future  vehicle  development.  Damage 
computations  will  also  be  conducted  to  demonstrate  how  the 
confined,  brittle  samples  behave.  The  study  shows  that  silicon 
carbide  and  boron  carbide  are  the  optimal  candidates  to  consider 
when  selecting  the  best  armor  performance  from  the  four 
configurations.  A  numerical  comparison  was  made  between  a 
pyroceram  confined  and  unconfined  configuration  and  ascertains 
approximately  a  twelve  percent  increase  in  ballistic  performance  of 
the  confined  sample.  The  computations  will  offer  the  researcher 
data  to  accurately  formulate  armor  to  protect  the  survivability  of  the 
ground  vehicle,  and  most  importantly,  the  soldier. 


INTRODUCTION 

Development  of  lightweight,  high  performance  armor  for  ballistic  protection  is  a 
critical  need  for  soldier  and  ground  platform  protection.  Experimental  range  testing  of 
projectiles  and  ballistic  armor  can  be  expensive,  time-consuming,  and  an  environmental 
hazard;  however,  it  is  necessary  to  verify  the  actual  behavior  of  such  systems.  Modeling 
and  simulation  (M&S)  allows  the  researcher  to  analyze  design  conditions  in  a  highly 
interactive  and  visual  environment  and  gives  details  about  the  impact  process  which 
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would  be  difficult  to  experimentally  obtain.  M&S  also  offers  specifics  into  the 
penetration  process  and  damage  accumulation  and  provides  information  for  the 
development  of  an  optimized,  experimental  test  plan,  reducing  both  cost  and  schedule. 
This  computational  study  quantifies  the  ballistic  efficiency  for  four,  confined  brittle 
materials  and  gives  an  understanding  about  the  physical  process  taking  place.  This 
analysis  will  also  compare  an  unconfined  and  confined  sample  and  evaluate  and  compare 
the  response. 

Attractive  materials  for  use  in  armor  systems  are  ceramics  due  to  their  superior 
hardness  and  compressive  strength.  It  is  important  to  confine  the  ceramic  fragments  and 
to  absorb  the  kinetic  energy  of  the  projectile  and  the  ceramic  debris  during  target 
penetration  and  increase  the  ballistic  efficiency.  Confinement  is  also  useful  for  multiple 
shot  testing  and  recovery  work.  Ceramics  have  a  low  density,  which  in  turn  reduces  the 
weight,  high  bulk  and  shear  moduli  that  prevent  large  deformations,  and  a  high  yield  stress  to 
preserve  resistance  to  failure.  Unfortunately,  ceramics  are  brittle  materials  that  produce 
extensive  fragmentation  caused  by  tensile  waves  generated  by  the  compressive  waves 
reflected  from  the  free  surfaces.  Carrying  out  a  variety  of  parametric  studies,  in  which  only 
one  parameter  can  be  altered,  clearly  demonstrates  each  design  limitation  and  the  overall 
effects  of  the  process  of  penetration  mechanics  and  damage. 

MODELING  APPROACH 

The  Lagrangian  hydrocode,  EPIC,  was  use  to  perform  the  computations  in  this 
parametric  study.  Lagrangian  methods,  in  which  the  finite  difference  grid  is  fixed  with 
the  material,  are  more  frequently  utilized  for  this  particular  type  of  structural  evaluation. 
Free  and  contact  surfaces  between  various  materials  remain  distinct  throughout  the 
calculation  and  the  mesh  deforms  with  the  evolutionary  flow.  The  four  cylindrical, 
confined  ceramic  specimens  investigated  were  aluminum  nitride,  boron  carbide,  9606 
pyroceram,  and  silicon  carbide.  The  confining  sleeve  was  modeled  with  RHA  and  had  a 
diameter  and  length  of  5  and  50  mm,  respectively.  The  front  and  rear  plug  materials  were 
both  RHA  with  a  diameter  and  length  of  25  and  10  mm.  The  RHA  sleeve  was  welded  to 
correlate  with  experimental  approaches.  The  tungsten  impactor  of  interest  was  modeled 
as  a  blunt,  long  rod  with  a  diameter  of  2  mm  and  length  of  80  mm.  Figure  1  illustrates 
the  modeled  geometry  of  the  projectile  used  in  the  analysis.  An  initial  computation 
involving  a  SiC  sample  with  a  specific  geometric  configuration  was  performed  based  on 
the  work  by  Holmquist  to  compare  computed  penetration  results  with  experimental 
findings  (Holmquist  [2002])  and  gain  insight  into  the  penetration  process  using  confined 
samples. 

The  problem  was  assumed  to  be  2-D  axisymmetric,  using  the  Johnson-Holmquist 
I  (Johnson  [1992]),  Johnson-Holmquist  II  (Johnson  [1994]),  and  Johnson-Holmquist- 
Beissel  (Johnson  [2003])  constitutive  models.  The  JH-1  model  did  not  allow  for 
continuing  damage;  whereas,  the  JH-2  constitutive  equation  includes  the  effect  of  damage 
and  residual  material  strength  (as  a  function  of  the  representative  variables)  and  the 
resulting  bulking  during  the  compressive  failure  of  a  ceramic  material.  The  JHB  equation 
is  more  recently  developed  and  better  describes  the  hydrostatic  pressure  including  a  phase 
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change  from  additional  experimental  data  for  various  ceramics.  In  this  analysis,  the  JH-1 
model  was  used  for  the  carbide  and  pyroceram,  the  JH-2  model  was  employed  for  the 
boron  carbide,  and  the  JHB  model  was  applied  for  the  aluminum  nitride  and  silicone 
carbide.  Tables  I,  II,  and  III  below  show  the  material  model  constants  used  for  the  JH-1, 
JH-2,  and  JHB  ceramic  cases. 
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Figure  1:  Original  configuration  (t  =  0  ps)  of  the  confined  samples 


Table  II.  Material  constants 
using  the  JH-2  constitutive 
model 

B4C 

Density  (Kg/m3) 

2510 

Specific  Heat  (J/kg-°C) 

999 

Shear  Modulus  (GPa) 

196 

Strength  Coefficient,  A 

0.93 

Strength  Coefficient,  N 

0.67 

Strength  Coefficient,  B 

0.70 

Strength  Coefficient,  M 

0.85 

Strain  Rate  Coefficient,  C 

0.005 

Tensile  Strength  (GPa) 

0.26 

Normalized  Fracture 
Strength  (GPa) 

0.20 

HEL  (GPa) 

19.0 

HEL  Pressure  (GPa) 

8.71 

HEL  Volumetric  Strain 
(GPa) 

0.041 

HEL  Strength  (GPa) 

15.4 

Damage  Constant,  D1 

0.001 

Damage  Constant,  D2 

0.50 

Bulk  Modulus,  K1  (GPa) 

233 

Pressure  Constant,  K2 
(GPa) 

-593 

Pressure  Constant,  K3 

(GP*) 

2800 

Table  III.  Material  constants 
using  the  JHB  constitutive  model 

AIN 

SiC 

Density  (Kg/m3) 

3226 

3215 

Specific  Heat  (J/kg-°C) 

735 

839 

Shear  Modulus  (GPa) 

127 

193 

Intact  Strength  (GPa) 

4.31 

4.92 

Pressure  (GPa) 

1.50 

1.50 

Maximum  Strength  (GPa) 

5.50 

12.2 

Strain  Rate  Coefficient,  C 

0.013 

0.009 

Failed  Strength  (GPa) 

0.09 

0.09 

Tensile  Strength,  T  (GPa) 

0.45 

0.75 

Maximum  Failed  Strength 
(GPa) 

0.20 

0.20 

Bulk  Modulus,  K1  (GPa) 

201 

220 

Pressure  Constant,  K2 
(GPa) 

260 

361 

Pressure  Constant,  K3 
(GPa) 

0 

0 

Damage  Constant,  D1 

0.16 

0.16 

Damage  Exponent,  N 

1.00 

1.00 

Table  I.  Material  constants  using  the  JH- 
1  constitutive  model 

Pyroceram 

Density  (Kg/m3) 

2608 

Specific  Heat  (J/kg-°C) 

975 

Shear  Modulus  (GPa) 

47 

Intact  Strength  1  (GPa) 

3.27 

Intact  Pressure  l(GPa) 

1.20 

Intact  Strength  2  (GPa) 

4.79 

Intact  Pressure  2  (GPa) 

4.12 

Strain  Rate  Coefficient,  C 

0.01 

Tensile  Strength  (GPa) 

0.24 

Maximum  Failed  Strength 
(GPa) 

3.11 

Bulk  Modulus,  K1  (GPa) 

76.0 

Pressure  Constant,  K2  (GPa) 

76.00 

Pressure  Constant,  K3  (GPa) 

0.000 

Damage  Constant,  D1 

4.1 

Damage  Constant,  EFMAX 

0.005 
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COMPUTED  RESULTS 

EPIC  was  accessed  and  the  analysis  was  carried  out  on  TARDEC’s  high 
performance  computing  (HPC)  system  on  a  32-processor  SGI®  Onyx®  350  for  the 
batch-oriented  computations.  The  HPC’s  hardware  was  used  to  provide  the 
computational  power  necessary  to  perform  this  analysis  and  save  time.  Figure  2  displays 
the  final  geometry  associated  with  the  impact  at  100  ps  into  the  four,  brittle  targets. 
EPIC  is  able  to  automatically  convert  the  2D  axisymmetric  finite  elements  into  meshless 
particles  as  they  become  severely  distorted,  which  can  be  easily  identified  in  the  figure. 
Majority  of  the  grid  remains  as  finite  elements,  and  this  allows  for  an  accurate  and 
efficient  solution  for  the  less  distorted  portion  of  the  problem  (Johnson  [2001]). 


AIN  B4C  Pyroceram  SiC 


Figure  2:  Final  geometry  associated  the  impact  velocity  of  1645  m/sec  into  the  four 
confined  specimens 

At  10  ps,  the  penetrator  completely  shears  through  the  RHA  front  plug  and 
reaches  the  face  of  the  ceramic  samples.  All  the  ceramics  begin  to  absorb  the  kinetic 
energy  and  fracture  and  lateral  cracking  of  the  ceramic  becomes  visible,  apart  from  the 
SiC  sample.  A  cavity  forms  in  the  targets,  the  long  rod  starts  to  erode  away,  and  the  spall 
plane  becomes  visible  at  the  front  surface.  The  ceramics  are  relatively  incompressible 
and  there  was  no  free  volume  for  the  comminuted  material  to  expand  into  except  the 
penetration  hole  itself.  The  rod  sheared  through  RHA  top  and  back  plate  by  90  ps  and 
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completely  exits  the  specimen  within  the  confined  perimeter  for  all  the  cases,  with  the 
exception  of  the  boron  carbide  and  silicon  carbide.  The  rod  completely  eroded  within 
and  is  arrested  by  the  B4C  target  configuration.  Interface  defeat  occurs  in  the  SiC 
sample,  or  the  condition  in  which  the  projectile  completely  dwells  on  the  surface,  but 
does  not  penetrate. 

Fig.  3  displays  the  damage  results  of  the  normal  impact  of  the  projectile  at 
50  ps  into  the  confined  ceramic  targets  with  a  50  mm.  thick  ceramic  sample  at  1645 
m/sec. 


:on TOURS 

■  0.00E-0  0 

■  2. DDE -01 

■  4.00E-01 


S.00E-01 
1  .OOE-OO 


AIN 


B4C  Pyroceram 


SiC 


Figure  3:  Damage  accumulated  from  the  impact  velocity  of  1645  m/sec  into  the  four 
confined  samples 


The  magnitude  of  the  compressive  stress  wave  exceeds  the  local  dynamic  strength 
of  the  material,  and  damage  begins  to  accumulate  through  the  formation  of  cracks  under 
the  impact  location.  By  10  ps,  a  comminuted  zone  is  formed  ahead  of  the  tungsten 
projectile  rod  due  to  the  intense  stress  conditions.  The  fracture  front  travels  through  the 
material  and  a  visible  conoid  of  pulverized  material  begins  to  form  under  the  impact 
location  for  the  first  four  samples,  with  the  exception  of  the  SiC  target.  It  is  obvious  the 
silicon  carbide  sample  has  the  best  ballistic  efficiency.  A  second  or  third  choice  material 
for  armor  solutions  would  be  boron  carbide  or  aluminium  nitride,  respectively. 

A  time  history  ASCII  file  was  updated  to  obtain  the  data  corresponding  to  the  z- 
coordinate  of  the  penetrator  in  order  to  trace  the  nose  into  the  target.  Figure  4  shows  the 
time  history  of  the  data  acquired  from  the  computations  performed. 
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100 


Figure  4:  Time  history  plot  of  the  nose  coordinates  for  the  long  rod  into  the  targets 
for  the  four  confined  configurations 

It  is  important  to  achieve  optimum  ballistic  performance  of  ceramic  targets  to  aid 
in  the  evaluation  of  the  response.  The  pyroceram  configuration  was  chosen  to  compare 
confinement  and  unconfined  targets,  because  the  tungsten  penetrator  remains  mostly  rigid 
during  the  computation.  Table  IV  shows  the  tungsten  rod’s  nose  coordinates  into  the 
actual  ceramics.  Figures  5  (a)  and  (b)  below  demonstrate  the  penetration  effects  and 
damage  accumulation  of  a  confined  and  unconfined  specimen  from  10  to  50  ps, 
respectively.  It  is  obvious  that  confining  the  ceramic  significantly  increases  the  ballistic 
efficiency.  By  60  ps,  the  rod  completely  exits  the  rear  of  the  unconfined  sample; 
whereas,  the  1645  m/s  threat  is  completely  arrested  in  the  confined  ceramic.  It  is  obvious 
that  confining  configurations  improve  the  mechanical  properties  by  inhibiting  the  failure 
mechanisms.  The  restraints  limit  the  flow  of  pulverized  ceramic  debris  and  can  be  used 
to  increase  the  resistance  to  penetration.  There  is  approximately  a  12%  increase  in 
ballistic  efficiency  for  the  confined  sample. _ 


Table  IV:  Nose  coordinates  of 
tungsten  rod  into  actual  ceramics 

Time  (|rS) 

Confined 

sample 

(mm) 

Unconfined 

sample 

(mm) 

0 

0 

0 

10 

7.772 

9.552 

20 

16.489 

18.389 

30 

24.593 

27.518 

40 

33.139 

36.579 

50 

41.854 

45.794 
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Figure  5:  Damage  progression  in  (a)  confined  and  (b)  unconfined  pyroceram 
computation  at  1645  m/s 
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SUMMARY  AND  CONCLUSION 


Computations  of  a  tungsten  long-rod  penetrator  impacting  four  confined  ceramics 
at  high-velocity  impact  were  conducted  to  propose  a  baseline  to  compare  ballistic 
performance.  Modeling  and  simulation  of  ballistic  events  offers  an  advanced 
understanding  of  the  failure  mechanisms  and  is  important  for  developing  improved 
models  and  for  designing  better  armor  systems.  Numerical  simulation  also  makes  the 
requirement  to  produce  hardware  for  experimental  testing  and  validation  more  efficient. 
In  this  study,  EPIC  was  employed  to  capture  the  final  geometry  and  damage  associated 
with  the  four  confined,  brittle  samples.  Confined  silicon  carbide,  boron  carbide,  or 
aluminum  nitride  are  the  optimal  candidates  to  consider  when  selecting  the  best  armor 
performance.  The  response  of  a  confined  and  unconfined  pyroceram  sample  was 
compared,  and  the  results  show  there  is  a  significant  increase  in  the  ballistic  performance 
of  the  confined  ceramic.  Additional  study  of  such  confined  configurations  is  needed  and 
could  provide  valuable  information  in  guiding  the  development  of  vehicle  armor  designs. 


REFERENCES 

Johnson,  G.R.  and  Holmquist,  T.J.  A  Computational  Constitutive  Model  for  Brittle 
Material  Subjected  to  Large  Strains,  High  Strain  Rates,  and  High  Pressures.  Shock- 
wave  and  High  Strain-rate  Phenomena  in  Materials,  1992. 

Johnson,  G.R.  and  Holmquist,  T.J.  An  Improved  Computational  Constitutive 
Model  for  Brittle  Materials.  High-Pressure  Science  and  Technology,  1993.  American 
Institute  of  Physics  Press,  1994. 

Johnson,  G.R.,  Stryk,  R.A,  Beissel,  S.R.,  and  Holmquist,  T.J.  Conversion  of 
Finite  Elements  into  Meshless  Particles  for  Penetration  Computations  Involving 
Ceramic  Targets.  Shock  Compression  of  Condensed  Matter,  2001. 

Holmquist,  T.J.  and  Johnson,  G.R.  Response  of  Silicon  Carbide  to  High  Velocity  Impact. 
Journal  of  Applied  Physics,  91,  2002. 

Johnson,  G.R.,  Holmquist,  T.J.,  and  Beissel,  S.R.  Response  of  Aluminum  Nitride 
(Including  a  Phase  Change)  to  Large  Strains,  High  Strain  Rates  and  High  Pressure. 
Journal  of  Applied  Physics,  99,  2003. 

Johnson,  G.R.,  Beissel,  S.R.,  Gerlach,  C.A.,  Stryk,  R.A.,  Holmquist,  T.J.,  Johnson,  A.A., 
Ray,  S.E.,  and  Arata,  J.J.  User  Instructions  for  the  2006  Version  of  the  EPIC  Code, 
Network  Computing  Services,  Inc.,  Minneapolis,  Minnesota,  2006. 


8 


